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We present a general derivation of semi-fermionic representation for spin operators in 
terms of a bilinear combination of fermions in real and imaginary time formalisms. The 
constraint on fermionic occupation numbers is fulfilled by means of imaginary Lagrange 
multipliers resulting in special shape of quasiparticle distribution functions. We show 
how Schwinger-Keldysh technique for spin operators is constructed with the help of 
semi-fermions. We demonstrate how the idea of semi-fermionic representation might be 
extended to the groups possessing dynamic symmetries (e.g. singlet /triplet transitions 
in quantum dots). We illustrate the application of semi-fermionic representations for 
various problems of strongly correlated and mesoscopic physics. 

PACS numbers: 71.27.+a, 75.20.Hr 

Introduction 

It is known that spin operators satisfy neither Fermi nor Bose commutation rela- 
tions. For example, the Pauli matrices for 5 = 1/2 operator commute on different 
sites and anticommute on the same site. The commutation relations for spins are 
determined by SU (2) algebra, leading to the absence of a Wick theorem for the gen- 
erators. To avoid this difficulty and construct a diagrammatic technique and path 
integral representation for spin systems various approaches have been used. The 
first class of approaches is based on representation of spins as bilinear combination 
of Fermi or Bose operators 1 - 7 , whereas the representations belonging to the second 
class deal with more complex objects like, e.g. the Hubbard 8 and supersymmetric 9 
operators, the nonlinear sigma model 10 etc. However, in all cases the fundamental 
problem which is at the heart of the difficulty is the local constraint problem. To 
illustrate it, let's consider e.g., first class of representations. Introducing the aux- 
iliary Fermi or Bose fields makes the dimensionality of the Hilbert space, where 
these operators act, greater than the dimensionality of the Hilbert space for the 
spin operators. As a result, the spurious unphysical states should be excluded from 
the consideration which leads in turn to some restrictions (constraints) on bilin- 
ear combinations of Fermi/Bose operators, resulting in substantial complication of 
corresponding rules of the diagrammatic technique. The representations from the 
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second class suffer from the same kind of problem, transformed either into a high 
nonlinearity of resulting model (non-linear sigma model) or hierarchical structure 
of perturbation series in the absence of Wick theorem (Hubbard operators). The 
exclusion of double occupied and empty states for a 5 = 1/2 impurity interact- 
ing with conduction electron bath (single impurity Kondo model), is controlled by 
fictitious chemical potential (Lagrange multiplier) of Abrikosov pseudofermions 5 . 
At the end of calculations this "chemical potential" A should be put A — > — oo to 
"freeze out" all unphysical states. In other words, there exists an additional U(l) 
gauge field which freezes the charge fluctuations associated with this representa- 
tion. The method works for dilute systems where all the spins can be considered 
independently. Unfortunately, attempts to generalize this technique to the lattice of 
spins results in the replacement of the local constraint (the number of particles on 
each site is fixed) by the so-called global constraint where the number of particles 
is fixed only on an average for the whole crystal. There is no reason to believe that 
such an approximation is a good starting point for the description of the strongly 
correlated systems. Another possibility to treat the local constraint rigorously is 
based on Majorana fermion representation. In this case fermions are "real" and cor- 
responding gauge symmetry is Z-i- The difficulty with this representation is mostly 
related to the physical regularization of the fluctuations associated with the discrete 
symmetry group. 

An alternative approach for spin Hamiltonians, free from local constraint prob- 
lem, has been proposed in the pioneering paper of Popov and Fedotov 11 . Based 
on the exact fermionic representation for 5=1/2 and 5=1 operators, where 
the constraint is controlled by purely imaginary Lagrange multipliers, these au- 
thors demonstrated the power and simplification of the corresponding Matsubara 
diagram technique. The semi-fermionic representation (we discuss the meaning of 
this definition in the course of our paper) used by Popov and Fedotov is neither 
fermionic, nor bosonic, but reflects the fundamental Pauli nature of spins. The goal 
of this paper is to give a brief introduction to a semi-fermionic (SF) approach. A 
reader can find many useful technical details, discussion of mathematical aspects of 
semi-fermionic representation and its application to various problems in the orig- 
inal papers n - 24 . However, we reproduce the key steps of important derivations 
contained in 19 , 20 in order to make the reader's job easier. 

The manuscript is organized as follows: in Section I, the general concept of 
semi-fermions is introduced. We begin with the construction of the SF formalism 
for the fully antisymmetric representation of SU (N) group and the fully symmetric 
SF representation of SU(2) group using the imaginary-time (Matsubara) represen- 
tation. We show a "bridge" between different representations using the simplest 
example of 5 = 1 in SU(2) and discuss the SF approach for 50(4) group. Finally, 
we show how to work with semi-fermions in real-time formalism and construct the 
Schwinger-Keldysh technique for SF. In this section, we will mostly follow original 
papers by the author 12 , 19 . The reader acquainted with semi-fermionic technique 
can easily skip this section. In Section II, we illustrate the applications of SF formal- 
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ism for various problems of condensed matter physics, such as ferromagnetic (FM), 
antiferromagnetic (AFM) and resonance valence bond (RVB) instabilities in the 
Heisenberg model, Dicke model, large negative - U Hubbard and t-J models, com- 
petition between local and non-local correlations in Kondo lattices in the vicinity 
of magnetic and spin glass critical points, dynamical symmetries in quantum dots, 
spin chains and ladders. In the Epilogue, we discuss some open questions and per- 
spectives. Necessary information about dynamical groups is contained in Appendix 
A, the effective models for spin chains are discussed in Appendix B. 

1. Semi-fermionic representation 

To begin with, we briefly reproduce the arguments contained in the original paper of 
Popov and Fedotov. Let's assume first S = 1/2. We denote as H the Hamiltonian of 
spin system. The standard Pauli matrices can be represented as bilinear combination 
of Fermi operators as follows: 

a) - a]aj - b)bj, af^2a]b j , aj^2b) aj . (1) 

on each site i of the lattice. The partition function of the spin problem Z a is given 
by 

Z a = Trexp(-/3H a ) = i N Tr cxp{- (3 (H F + iirN F /(2(3)) (2) 
where Hp is the operator obtained from H a by the replacement (1) and 

N 

N = £(at 0j + b%) (3) 
j=i 

(N is the number of sites in the system and j3 = 1/T is inverse temperature). To 
prove equation (2) we note that the trace over the nonphysical states of the i-th 
site vanishes 

Tr unphys exp(-/3(ff F + inN F /(20)) = (-*) + {-if = (4) 

Thus, the identity (2) holds. The constraint of fixed number of fermions Nj = 1, 
is achieved by means of the purely imaginary Lagrange multipliers \i = —itt/{2(3) 
playing the role of imaginary chemical potentials of fermions. As a result, the Green's 
function 

G=(iuj F -e)- 1 (5) 

is expressed in terms of Matsubara frequencies ui F — 2nT(n +1/4) corresponding 
neither Fermi nor Bose statistics. 

For S — 1 we adopt the representation of H a in terms of the 3-component Fermi 
field: 

u] -> a]a 3 - b]b jt a+ -» V2(a] Cj + c}b 3 ), aj - V2(c]a 3 + b\ Cj ). (6) 
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The partition function Z a is given by 



Z a = Tr(-f3H a ) = ( -L ) Tr exp(-0(H F + iirN F /(3f3)). 




(7) 



It is easy to note that the states with occupation numbers and 3 cancel each 
other, whereas states with occupation 1 and 2 are equivalent due to the particle- 
hole symmetry and thus can be taken into account on an equal footing by proper 
normalization of the partition function. As a result, the Green's function in the 
imaginary time representation is expressed in terms of lof = 2nT(n + 1/3) frequen- 
cies. 

In this section, we show how semi-fcrmionic (Popov-Fcdotov) representation 
can be derived using the mapping of partition function of the spin problem onto the 
corresponding partition function of the fermionic problem. The cases of arbitrary 
N (even) for SU(N) groups and arbitrary S for SU(2) group are discussed. 

1.1. SU(N) group 

We begin with the derivation of SF representation for SU(N) group. The SU(AT) al- 
gebra is determined by the generators obeying the following commutation relations: 



where a, (3 = 1,...,N. We adopt the definition of the Cartan algebra 25 of the SU(iV) 
group {H a } = S% similar to the one used in 26 , noting that the diagonal generators 
are not traceless. To ensure a vanishing trace, the diagonal generators should 
only appear in combinations 



which effectively reduce the number of independent diagonal generators to N — 1 
and the total number of SU(Y) generators to TV 2 — 1. 

In this paper we discuss the representations of SU(N) group determined by 
rectangular Young Tableau (YT) (see 26 and 19 for details) and mostly concentrate 
on two important cases of the fully asymmetric (one column) YT and the fully 
symmetric (one row) YT (see Fig.l). 

The generator Sp may be written as biquadratic form in terms of the Fermi- 




(8) 



N 



N 




(9) 



operators 




(10) 



7 



where the "color" index 7 = 1, ...,n c and the n c (n c + l)/2 constraints 



N 




(11) 
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m 



a) 



b) 



c) 



Fig. 1. Rectangular Young tableau to denote a SU(iV) representation (a), single column tableau 
corresponding to n c = 1 (b) and single row tableau standing for spin S = n c /2 representation of 
SU(2) group (c). 



restrict the Hilbert space to the states with m * n c particles and ensure the charac- 
teristic symmetry in the color index a. Here m corresponds to the number of rows 
in rectangular Young Tableau whereas n c stands for the number of columns. The 
antisymmetric behavior with respect to a is a direct consequence of the fermionic 
representation. 

Let us consider the partition function for the Hamiltonian, expressed in terms 
of SU(./V) generators 



Z s = Trcxp(-(3H S ) = Tr' exp(-/J# F ) 



(12) 



where Tr' denotes the trace taken with constraints (11). As it is shown in 19 , the 
partition function of SU (N) model is related to partition function of corresponding 
fcrmion model through the following equation: 



Z s = J IT d»U)P(Kj))Tr exp (-0(H F - M (j>f 



)) 



(13) 



here P{nj) is a distribution function of imaginary Lagrange multipliers. We calculate 
P([ij) explicitly using constraints (11). 

We use the path integral representation of the partition function 

Zs/Z% = J J] d»(j)P(ii(j)) cxp(A)/ J [] dn(j)P(nti)) cxp(A) (14) 

3 3 

where the actions A and Aq are determined by 

,/3 n p 

A = A Q - dTH F (r), A = VV dTa k {j,T){d T + fJ,(j))a k (j,T) (15) 

- k=1 Jo 
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and the fermionic representation of SU(iV) generators (10) is applied. 

Let us first consider the case n c = 1. We denote the corresponding distribution 
by PN,m(n(j))> where m is the number of particles in the SU(A r ) orbital, or in other 
words, 1 < m < N labels the different fundamental representations of SV(N). 



N 



^2a k (j)a k (j) = m (16) 



fe=i 

To satisfy this requirement, the minimal set of chemical potentials and the corre- 
sponding form of PN, m {n{j)) are to be derived (see Fig. 2). 

To derive the distribution function, we use the following identity for the con- 
straint (16) expressed in terms of Grassmann variables 

s- 1 • / l w / • ( ^{ n i - m )\ 

<^, m = — sin(7r(nj-m))/sinl J — I (17) 

Substituting this identity into (12) and comparing with (14) one gets 

PnMD) = x E CX P HvT^ !) ) S <J*U) Wt), (18) 
fe=i ^ ' 

where 

V k = - — (2k-l). (19) 

Since the Hamiltonian is symmetric under the exchange of particles and holes when 
the sign of the Lagrange multiplier is also changed simultaneously, we can simplify 
(18) to 

2 i L JV / 2 J / 2 fc-l\ 
PN, m (p(j)) = jj E sin Urn— — J 5(p{j) - fi k ) (20) 
fe=i ^ ' 

where 2 J denotes the integer part of N/2. As shown below, this is the minimal 
representation of the distribution function corresponding to the minimal set of the 
discrete imaginary Lagrange multipliers. Another distributions function different 
from (20) can be constructed when the sum is taken from k = N/2 + 1 to N. 
Nevertheless, this DF is different from (20) only by the sign of imaginary Lagrange 
multipliers fi k = /i£ = — [i k and thus is supplementary to (20). 

Particularly interesting for even N is the case when the SU(iV) orbital is half- 
filled, m = N/2. Then all Lagrange multipliers carry equal weight 

2 . N / 2 

P N ,N/2(Kj)) = TfEC-i^W) - w0- (2i) 

fe=l 

Taking the limit N — > oo one may replace the summation in expression (21) in a 
suitable way by integration. Note, that while taking N — > oo and m — > oo lim- 
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Fig. 2. Graphical representation of exclusion principle for SU(iV) semi-fermionic representation 
with even N, n c = l (we use fj, = inT/2 for SU(2) and fn = inT/4, fi 2 = 3ivrT/4 for SU(4)). 



its, we nevertheless keep the ratio m/N = 1/2 hxed. Then, the following limiting 
distribution function can be obtained: 







N 



(22) 



resulting in the usual continuous representation of the local constraint for the sim- 
plest case n c = 1 



Z s = Tr(exp {-/3H F ) S ( rij 



N 
~2 



(23) 



We note the obvious similarity of the limiting DF (22) with the Gibbs canonical 
distribution provided that the Wick rotation from the imaginary axis of the Lagrange 
multipliers \i to the real axis of energies E is performed and thus [i(j)N/2 has a 
meaning of energy. 

Up to now, the representation we discussed was purely fermionic and expressed in 
terms of usual Grassmann variables when the path integral formalism is applied. The 
only difference from slave fermionic approach is that imaginary Lagrange multipliers 
are introduced to fulfill the constraint. Nevertheless, by making the replacement 

f iirr 2k — 1\ .. . . / MTT2fc— 1\. . 
ak\3,T)) -> a fe (j,r)exp — — - , a fe (j,r) -> a fc (j,r)exp — - (24) 



(3 N 



(3 N 



I 



l 
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we arrive at the generalized Grassmann (semi-fermionic) boundary conditions 



Ofc0',/3) = a k (j, 0)exp m 



2k - 1 



dk{j,P) = a k (j, 0) exp -Z7T 



2k- 1 



(25) 



This leads to a temperature diagram technique for the Green's functions 

G a(3 (j,r) = -(T T a a (j,T)a f3 (j,0)) (26) 

of semi-fermions with Matsubara frequencies different from both Fermi and Bosc 
representations (see Fig. 3). 




Fermi 



n + ^f n + K 1 -^) 




n+1/2 



2jiT 



Fermi 




n+ ; 



..n + |(l- 



-+ 
n+1/2 



C0„ 
2jiT 



Fig. 3. The minimal set of Matsubara frequencies for a) SU(N) representation with even N/ 
SU{2) representation for half-integer value of the spin, b) 5(7(2) representation for integer values 
of the spin and 1 = 1. 



The exclusion principle for this case is illustrated on Fig. 2, where the 5=1/2 
representation for the first two groups SU(2) and SU(4) are shown. The first point 
to observe is that the spin Hamiltonian does not distinguish the n particle and the 
n hole (or N — n particle) subspace. Eq. (19) shows that the two phase factors 
exp(/3/in) and exp(/3/i(iV — n)) accompanying these subspaces in Eq. (20) add up 
to a purely imaginary value within the same Lagrange multiplier, and the empty 
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and the fully occupied states are always cancelled. In the case of N > 4, where we 
have multiple Lagrange multipliers, the distribution function P{p) linearly combines 
these imaginary prefactors to select out the desired physical subspace with particle 
number n = m. 

In Fig. 2, we note that on each picture, the empty and fully occupied states are 
cancelled in their own unit circle. For SU(2) there is a unique chemical potential 
/i = ±iwT/2 which results in the survival of single occupied states. For SU(4) there 
are two chemical potentials (see also Fig. 3). The cancellation of single and triple 
occupied states is achieved with the help of proper weights for these states in the 
distribution function whereas the states with the occupation number 2 are doubled 
according to the expression (21). In general, for SU(iV) group with n c = 1 there 
exists N/2 circles providing the realization of the exclusion principle. 

1.2. SU(2) group 

We consider now the generalization of the SU(2) algebra for the case of spin S. 
Here, the most convenient fermionic representation is constructed with the help of 
a 2S + 1 component Fermi field a,k(j) provided that the generators of SU(2) satisfy 
the following equations: 

S-i 

5+ = ]T VS(S+l)-k(k + l)al +1 (j)a k (j), 
k=-s 

s 

S-= J2 A/5(S + l)-fe(fe-l)4-iC7'K(j), 

fe=-S+l 

s 

S z = £ ka\{j)a k {j) (27) 

fe=-S 

such that dimHp = 2 2S+1 whereas the constraint reads as follows 

k=S 

nj= E 4(JK0') = * = 1 (28) 

k=-S 

Following the same routine as for SU(AT) generators and using the occupancy con- 
dition to have I = 1 (or 25*) states of the (25*+ 1) states filled, one gets the following 
distribution function, after using the particle-hole symmetry of the Hamiltonian 

2l Ls i^ /2J / ik -w 

%+i,iW))=257i E s ™(*MTl) 6 (»ti)-^) (29) 

fe=i ^ ' 

where the Lagrange multipliers are Uk = — «7rT(2fc— 1)/(2S+ 1) and k = 1, [S + 
1/2J, similarly to Eq.(19). 
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In the particular case of the SU(2) model for some chosen values of spin 5 the 
distribution functions are given by the following expressions 



W)) = ««Ut?) + 



iwT\ 



for 5 = 1/2 



p 3 ,i(Mi)) = PsMo)) = j= s (Mi) + —) 



(30) 



(31) 



for 5 = 1. 

This result corresponds to the original Popov-Fedotov description restricted to 
the 5 = 1/2 and 5=1 cases. 




Fig. 4. Graphical representation of exclusion principle for SU(2) semi-fermionic representation 
for 5 = 1 and 5 = 2. For any arbitrary integer value of spin there exists S circle diagrams 
corresponding to the 5 different chemical potentials and providing the realization of the exclusion 
principle. 



We present as an example some other distribution functions obtained according 
to general scheme considered above: 



iV2 



6(n 



inT x 



5(fx 



3z7rT, 



(32) 
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for 5 = 3/2, SU(2) and 

PaAv)= 2 { 5 ^+— )-^+— )J ( 33 ) 
for effective spin "5 = 1/2", SU(4), 

^S,l(/i) = ^5,4 M = 

for 5 = 2, 5C/(2) etc. 

A limiting distribution function corresponding to Eq. (22) for the constraint 
condition with arbitrary I is found to be 

Poo,l(j*U)) ^ ^~ cxp(-^O')). (35) 

For the case I = m = N/2 — > oo and 5 = (N — l)/2 — > oo the expression for 
the limiting DF Poo.i(p(j)) coincides with (23). We note that in S — > oo (or N — > 
oo) limit, the continuum "chemical potentials" play the role of additional U(l) 
fluctuating field whereas for finite S and N they are characterized by fixed and 
discrete values. 

When S assumes integer values, the minimal fundamental set of Matsubara 
frequencies is given by the table in Fig. 3. 

The exclusion principle for SU(2) in the large spin limit can be also understood 
with the help of Fig. 2 and Fig. 4. One can see that the empty and the fully occupied 
states are cancelled in each given circle similarly to even-iV SU(iV) algebra. The 
particle-hole (PH) symmetry of the representation results in an equivalence of single 
occupied and 2S occupied states whereas all the other states are cancelled due to 
proper weights in the distribution function (29). In accordance with PH symmetry 
being preserved for each value of the chemical potential all circle diagrams (see Fig. 3, 
Fig. 5) are invariant with respect to simultaneous change /i <-» — fj, and n partic i e <-» 

Wholes • 

1.3. Semi-fermionic representation of SO(N) group 

In this chapter we show the way of generalization of semi-fermionic representation 
to the dynamical algebras o(4) and o(5) (see Appendix A). Like in the case of pure 
spin operators, these representations should preserve all kinematical constraints. 

The first step to derive semi-fermionic representation for SO (4) group is based 
on the local isomorphism of 50(4) and SU(2) x SU(2). 

We start with n = 4 - field representation of 5/7(2) group (10) 



(an, ai2, a2i, a 2 2) 



(36) 
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There are two diagonal and two off-diagonal constraints (11) which read as follows: 

a\l a lt + a 21 a 21 = 1, «l 2 a 12 + «22 a 22 = 1, (37) 

Onfli2 + a 2 i°22 = 0, a\ 2 an + a\ 2 a 2 i = 0, (38) 

and generators ol SU (2) group are given by 

S~ = S 2 = a} 1 a2i + a\ 2 a-22, S + — Sf — a' 21 au + a\ 2 a 12 

2S Z = S 2 — S\ = a\ x a 2 i + a\ 2 a 22 - oj^n - a\ 2 a 12 (39) 

Combining definition (39) with constraint (38) we reach the following equations: 

S~ = 0^(021 + 012) + (012 + 021)022, 
S+ = (021 + a i2) a n + a\ 2 (ai 2 + o 2 i), 

5 Z = a^ 2 a 22 — a^an (40) 

Therefore, we conclude that the antisymmetric (singlet) combination 012 — 021 does 
not enter the expression for spin 5 = 1 operators. Thus, three (out of four) com- 
ponent Fermi-field is sufficient for the description of S = 1 SU (2) representation in 
agreement with (27). Defining new fields as follows 

011 = f-i, a 22 = /1, 7^( ai2 + fl21 ) = f°> tJ| (° 12 - fl21 ) = A- ( 41 ) 

where fermions /i,/o,/-i stand for 5 Z = 1,0 — 1 projections of the triplet state 
and fermion s determines the singlet state, we come to standard S = 1 SU{2) 
representation 

5+ = >/2(/ V-i + AVo), = ^(/li/o + /0V1), 5 Z = AVi - /li/-i, (42) 

with the constraint 

ni + n + n_i + ra s = 2 (43) 

where n Q = /t/ a . 

The constraint (43) transforms to a standard SU(2) S — 1 constraint (48) in 
both cases n s = and n s = 1 since there is no singlet/triplet mixing allowed by 
SU(2) algebra. 

To demonstrate the transformation of the local constraint let's first consider the 
case n s = 0. The constraint reads as follows 

ni + n + n_i = 25 <^=> S 2 = 5(5 + 1). (44) 

On the other hand, the states with 25 occupation are equivalent to the states with 
single occupation due to particle-hole symmetry. Thus, the constraint (43) might 
be written as 

n\ + h + n_i = 1 (45) 
where n a = 1 — n a . The latter case corresponds to n s = 1. 
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The singlet /triplet mixing is allowed for 50(4) group. This mixing is described 
in terms of 3 additional generators responsible for transitions between singlet and 
triplet (see Appendix A). Using 4-component auxiliary fermions fx, where A = 
— 1,0, l,s and matrix form of generators (A. 13) we represent R - generators of 
50(4) as follows 

R+ = V2(ftf s - /]/_!), R~ = x/2(/j/i - /li/s), R z = -(/oV. + /]/o). 

with the only constraint 

m +n +n_i +n s = 1, ^> S 2 + R 2 = 3 (46) 

whereas the orthogonality condition S-R = R- S = 0is fulfilled automatically. The 
constraint (46) is respected by means of introducing real chemical potential A — > oo 
for Abrikosov's auxiliary fermions or imaginary chemical potentials fit — —inT/S 
for Popov-Fedotov semi-fermion. 

The fermionic representation of 50(5) group is easily constructed by use semi- 
fermionic representation and is characterized by 5- vector q T = (/Ij/q, fl,fj,fr) 

5+ = v^/oV-i + flfo), S* = flh - fUf-i, 

R+ = V2(ftf s - /]/_!), R* = -(/t/ s + /t/ ), 

P+ = V2(.flf r - /+/_!), P z = -(flu + /t/ ). (47) 

and 

A = i{ftfs flfr) 

The constraint 

ni + n + n_i + n s + n r = 1 (48) 

is respected either by real infinite chemical potential (Abrikosov pseudofermions) 
or by set of complex chemical potentials (semi- fermions). We do not present here 
ten 5x5 matrices characterizing 50(5) representation to save a space. The reader 
can easily construct them using representations (47). There exists also a bosonic 
representation based on Schwinger bosons which might be derived by the method 
similar to used above for 50(4) group. The representations of higher 50(n) groups 
can be constructed in a similar fashion. 

Kinematic constraints imposed on auxiliary fermions and bosons is in strict com- 
pliance with the Casimir and orthogonality constraints in spin space. Accordingly, 
the number of fermionic and bosonic fields reproduces the dimensionality of spin 
space reduced by these constraints. We have seen that the 6-D space of generators of 
50(4) group is reduced to D—A. Then the minimal (unconstraint) fermionic repre- 
sentation for this group should contain two U(l) fermions. This means that the rep- 
resentation (46) is not minimal. Apparently, the best way to find such representation 
is to use Jordan- Wigner-like transformation 10 . The exact form of Jordan- Wigncr 
transformation for 50(4) group was recently reported in 27 . The same kind of ar- 
guments applied to 50(5) group tells us that the spinor field should contain seven 
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components. This means that the 3-color U(l) fermionic representation should be 
completed by one more real (Majorana) lermion, and this fact points to one more 
hidden Z 2 symmetry 10 . 



1.4. Real-time formalism 

We discuss finally the real-time formalism based on the semi-fcrmionic represen- 
tation of SU(N) generators. This approach is necessary for treating the systems 
out of equilibrium, especially for many component systems describing Fermi (Bose) 
quasiparticles interacting with spins. The real time formalism 28 , 29 provides an al- 
ternative approach for the analytical continuation method for equilibrium problems 
allowing direct calculations of correlators whose analytical properties as function of 
many complex arguments can be quite cumbersome. 

To derive the real-time formalism for SU(iV) generators we use the path inte- 
gral representation along the closed time Keldysh contour (see Fig. 5). Following 



OO 



t-i p 



D + OO 



Fig. 5. The Keldysh contour going from — oo — > oo — * — oo in real time. The boundary conditions 
on the imaginary time segment determine the generalized distribution functions for quasiparticles. 



the standard route 30 , we can express the partition function of the problem con- 
taining SU(iV) generators as a path integral over Grassmann variables tpi = 
(ais(j), ai,N(j)) T where I = 1, 2 stands for upper and lower parts of the Keldysh 
contour, respectively, 

Z/Z = J D^D^exp(iA)/ J D$Dipexp(iA ) (49) 

where the actions A and Ao are taken as an integral along the closed-time contour 
Ct + C T which is shown in Fig. 5. The contour is closed at t = — oo + it since 
cxp(— /3Hq) = T T exp ^— Jq Hq(1t^ . We denote the ip fields on upper and lower 
sides of the contour Ct as tpi and ip2 respectively. The fields ^> stand for the contour 
C T . These fields provide the matching conditions for i/ji j2 and are excluded from 
the final expressions. Taking into account the semi-fermionic boundary conditions 
for generalized Grassmann fields (25) one gets the matching conditions for tpi^ at 
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t = ±00, 

i>l, a \k(-00) = CXp (iK 2k N ~ ) ^cUC" 00 ), 

<aU(+°°)=1#,alfc(+°°) (50) 

for k = 1, L-^/2j and a = 1,...,N. The correlation functions can be represented 
as functional derivatives of the generating functional 

Zfa] = Z o 1 J Di>DiP exp (iA + ij> dt(fja z ^ + il>a z r])j (51) 

where 77 represents sources and the a z matrix stands for " causal" and " anti-causal" 
orderings along the contour. 

The on-site Green's functions (GF) which are matrices of size 2N x 2N with 
respect to both Keldysh (lower) and spin-color (upper) indices are given by 

^ M ~- i imidm m ^ (52) 

To distinguish between imaginary-time (26) and real-time (52) GF's, we use different 
notations for Green's functions in these representations. 

After a standard shift-transformation 30 of the fields ip the Keldysh GF of free 
scmi-fermions assumes the form 

rial ,\ _ /-tR,a ( 1 — fe ~fe \ /ut,o ( ~ fe ~ fe \ 

Go(e) - G ° U-/.-/J" U-/«i-/J' 

where the retarded and advanced GF's are 

G^' A)a (e) = (e±^)- 1 , / e = /W fc )(e), (53) 
with equilibrium distribution functions 

1 ^ iun k - e e^(2fe-i)/JV exp (/3e) + 1 ' 1 ' 

A straightforward calculation of f( N > k ) for the case of even N leads to the following 
expression 

N 



^(-l)'- 1 cxp(/3e(iV-0)exp(- 
i=i ^ 



M(2k- 1) 



N 

^"•"w = ^ (55 » 

where fc = I,..., N/2. The equilibrium distribution functions (EDF) f( 2S + 1 . k ) fo r 
the auxiliary Fermi-fields representing arbitrary 5 for SU(2) algebra are given by 

^(-l)'- 1 exp (/3e(2S + 1 - I)) exp (~ )) 
/ (2S+i,fc) (c) = _ exp((25+l)/3e) + (-l) 2S + 1 (56) 



I 



I 
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for k = 1, [S + 1/2J . Particularly simple are the cases of 5* = 1/2 and S = 1, 

"'<*>- 4 555(j»o <57) 

,<»><.)- -£..<*, -Wi*^ (58, 

Here, the standard notations for Fermi/Bose distribution functions n F / B (e) = 
[exp(/3e) ± 1] _1 are used. For 5=1/2 the semi-fermionic EDF satisfies the ob- 
vious identity |/ (2,1) (e)| 2 = n F (2e). 

In general the EDF for half-integer and integer spins can be expressed in terms 
of Fermi and Bose EDF respectively. We note that since auxiliary Fermi fields intro- 
duced for the representation of SU(iV) generators do not represent the true quasi- 
particles of the problem, helping only to treat properly the constraint condition, the 
distribution functions for these objects in general do not have to be real functions. 
Nevertheless, one can prove that the imaginary part of the EDF does not affect the 
physical correlators and can be eliminated by introducing an infinitcsimally small 
real part for the chemical potential. In spin problems, a uniform/staggered magnetic 
field usually plays the role of such real chemical potential for semi-fcrmions. 



2. Application of semi-fermionic representation for strongly 
correlated systems 

In this Section we illustrate some of the applications of SF representation for various 
problems of strongly correlated physics. 



2.1. Heisenberg model: FM, AFM and RVB 

The effective nonpolynomial action 31,32 for Heisenberg model with ferromagnetic 
(FM) coupling has been investigated in 11 . The model with antiferromagnetic 
(AFM) interaction has been considered by means of semi-fermionic representation in 
17 and 18 (magnon spectra) and in 12 for resonance valence bond (RVB) excitations. 
The Hamiltonian considered is given as 

H mt = - £ Jij Usj \) (59) 

<ij> ^ ' 

• Ferromagnetic coupling J = Ifm > 

The exchange SiSj is represented as four-semi-fcrmion interaction. Applying the 
Hubbard-Stratonovich transformation by the local vector field <J>i(r) the effective 
nonpolynomial action is obtained in terms of vector c-field. The FM phase transition 
corresponds to the appearance at T < T c of the nonzero average (<E> 2 (q = 0, 0)) which 
stands for the nonzero magnetization, or in other words, corresponds to the Bose 
condensation of the field $ z . 

$ z (k, u) = -M(/3^) 1/2 %,o<^o + &(k, w). (60) 
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In one loop approximation the standard molecular field equation can be reproduced 

M = I FM (0)tanh(/3M/2). (61) 
The saddle point (mean-field) effective action is given by well-known expression 

f3M 2 . / , ff3M' 



A [M] = -N 



4/ M (0) 



In 2 cosh 



and the free energy per spin f is determined by the standard equation: 



\nZ s 



pM 2 
4/ M (0) 



— In 2 cosh 



0M 



(62) 



(63) 



Calculation of the second variation of A e f / gives rise to the following expression 



2cosh 2 (/3rj) 
"7 E lM{k)&>(k,u)**(k,u)- 



$ z (fc,0)- 



k,u>=£0 



■E $+ &<") 



tanh(/?ft) 
2S1 — 



$ _ (A,w) 



(64) 



where = (gusH + A4)/2. The magnon spectrum (T < T c ) is determined by the 
poles of ($ + $~) correlator, w = Ak 2 . 

• Antifcrromagnetic coupling J = Iafm < 0. Ned solution 

The AFM transition corresponds to formation of the staggered condensate 

**(£,o;) = tf(pN) 1/ %,tf u ,o + *'$,<*>) (65) 

The one-loop approximation leads to standard mean-field equations for the stag- 
gered magnetization 



A [M] = N 



N= -/4f M (Q)tanh(/W/2), 
/W 2 



In 2 cosh 



m 



(66) 



_4/afm(<2) 

After taking into account the second variation of A e ff, the following expression 
for the effective action is obtained [(see e.g. 17 , 18 ): 



2cosri 

a; 



$ Z (M)+ 
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k.uj 



AFM 



(*) + 



2fhanh(/3f2) 
4ft 2 + u> 2 



$ _ (jfe,w)- 



The AFM magnon spectrum cj = c|k|. 



4ft 2 + 



(67) 



• Antifcrromagnetic coupling. Resonance Valence Bond solution 

The four-semi-fermion term in (59) is decoupled by bilocal scalar field Ay. The 
RVB spin liquid (SL) instability in 2D Heisenberg model corresponds to Bose- 
condensation of exciton-like pairs of semi-fermions (for simplicity we consider the 
uniform RVB state) : 



i0 = _ E i tMi (i£.), 

q 

Ao = - ln I 2 cosh(/5/ q A ) 



(68) 



where A = A(q = 0) is determined by the modulus of Ay field 



A 



<ij> 



(R, f) = A(r) exp [irA(R) 



(69) 



whereas the second variation of SA e ff describes the fluctuations of phase Ay 

A eff = £ A»(k,wKj,^/8(k,w), 



a/3 



Tr(pV(G p+fc G p + G p+fe G p ) + <5 Q/3 /(/ p A )) 



(70) 



The spectrum of excitation in uniform SL is determined by zeros of ir R and is 
purely diffusive 33 - 34 . The recent development of application of semi-fermions to 
low dimensional magnetic systems can be found in 21 > 22 . 



2.2. Dicke model 

In this Section we describe the application of semi-fermionic approach to two-level 
systems interacting with single-mode radiation field (Dicke model). The influence 
of dissipative environment on two-level system has been extensively studied in spin- 
boson model 35 (see also a review 36 for coherent effects in mesoscopic few level 
system, Dicke super and sub-radiance effects) . Addressing the reader to above men- 
tioned reviews for discussion of physical implementation of the Dicke model, we 
discuss in this section only technical aspects related to derivation of its equilibrium 
(thermodynamical) properties. We closely follow original derivation contained in 
Popov- Fedotov paper 11 . 
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The Dicke Hamiltonian 

O N N 
H D = too^ip + + £ + V>V) (71) 

i i 

contains cr-matrices representing two-level systems and bosonic V'-operators describ- 
ing the single-mode radiation field. 

Applying semi-fermionic representation for two-level system one gets 

N N 
H D = u^H + - ]T (ajoi - b\bi) + E + ^ b l a ) ( 72 ) 



where a, and bi stand for semi-fermionic fields with generalized Grassmann bound- 
ary conditions. 

The semi-fermionic variables appear quadratically in the action and can be inte- 
grated out. As a result, the partition function is represented as a ratio of two path 
integrals 

Z /Z „ = J D[^} C xp(So[mdctMm)) N / J Dty] exp(S [V])(det M([0})) N (73) 
where 

So W\ = E («w - wo) i>* M^(w) 

UJ 

and M is an operator with elements 

_/ (zp + n/2)5 m g/((3N)V*r(p-q)\ (7A] 
Mpq -{g/(f3N) 1/2 ^q-p) (ip-mS Pq )' { ] 

Here uj = 2irnT is bosonic and p — 2-nT(n +1/4) is semi-fermionic Matsubara 
frequencies. 

Evaluating the integrals (73) on gets the following asymptotic for partition func- 
tion: 

1 



+ 0(1/N), T>T C (75) 



1 — a(uj) 

Z Q /Z Qa =AN 1 ^ 2 eMBN)Y[L-\uj) + 0(l/N 1 / 2 ), T <T C (76) 

UJ 

where T c is determined from the equation 

' / ! tanh (SL\ = l (77) 



uj Q \2T c 
and the following notations are used: 



(UJo - lLU)(il — lUJ) 



g 2 (l- pn A /smh(l3n A )) 



i 



i 
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B = m[cosh(/3Q A /2)/cosh(/?ft A /2)] - uj A 2 P, n 2 A = Q 2 + 4g 2 A 2 . (80) 

The parameter A(T) is determined from Eq.77 with the replacement 17 — > f2 A and 
T C -»T. 

The Bose spectra are obtained by analytic continuation iuo — > E + iS of the 
equations 

1 - a(uj) = 0, T > T c L(lu) = 0, T < T c (81) 

and are given by 

Ei,2 = \ 



fi + woT - ^o) 2 + V tanh(/30/2) , T > T c (82) 

£i=0, S 2 = [[(r! + wo) 2 +4g 2 A 2 ] 1/2 , T < T c . (83) 

Multimode variants of the Dicke type and also Dicke models with interaction that 
takes non-resonance term into account can be investigated analogously. 

2.3. Spin- dependent semi-fermionic representation for Hubbard 
and t-J models 

In this Section we discuss the spin-dependent Popov-Fedotov representation and 
application of semi-fermions to Hubbard and t-J models. 

We start with the negative - U Hubbard model described by the Hamiltonian 

H = - ^ tijal a a iia + U^2 n A n i\- ( 84 ) 

The interaction U < is attractive and the hopping amplitude is assumed small 
\t\ <C \U\. The physical situation underlining this limit of model is characterized 
by the existence of exactly one single occupied site and empty or double occupied 
sites otherwise. The physical restriction also imposes a subtle constraint on the 
corresponding Hilbert space, which requires that all contributions from states with 
more than one unpaired electron are ruled out. Thus, the situation is opposite the 
Popov-Fedotov limit where only single occupied states represent physical states of 
the model, while empty and double occupied states were unphysical and eliminated 
by the single occupancy condition. 

In order to remove single occupied states in large negative-U Hubbard model 
(LNU) from thermodynamic averages, it was proposed 14 ' 37 to introduce an imagi- 
nary magnetic field to the Hamiltonian 

H = H -ianT^in^ -n n ) (85) 

i 

with a as a real parameter (—1 < a < 1). The purely imaginary magnetic field term 
does not violate time-reversal symmetry 14 . 

The local partition function = Tr; t cxp(—f3H) is given by 

ZW = (1, l|e-^|l, 1) + (0, 0|e-^|0, 0) + 7 (1, 0| e -^|l, 0) (86) 
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where symmetry condition 

(i,o|ff|i,o) - (o,i|#|o,i) 

is used and 7 = 2 cos(-7ra) is a statistical parameter. The lull lermionic Hilbcrt space 
is recovered if a = (7 = 2). The case a = ±1/2 (7 = 0) result in the semi-fermionic 
theory for the LNU-Hubbard model. 

If one employs the standard Matsubara temperature diagram technique, it can 
be easily seen that the additional term in the Hamiltonian H enters in the one- 
particle Green's function as a spin-dependent imaginary chemical potential, which 
alternatively may be absorbed in the Matsubara frequencies 

w„ = 2nT(n + [l + a<r]/2) 

with a — ±1. These spin-dependent frequencies interpolate between the bosonic 
and fcrmionic result and constitute the semi-fermionic description of the model. 

The mean- field result for the occupation probability (n icr ) wit the chemical po- 
tential (i = fi + U — Azt 2 /Uu (z is a coordination number, v is a filling factor) gives 
the complex value 

M = cxp(-/^) + 1 + *2cosh(/3 M /2) (8?) 

similarly to complex distribution function derived in real-time spin formalism (57). 
The result reflects the fact that he electron occupation probability (n^) is not an 
observable quantity in the strong coupling limit, whereas its real part Re(n^} — 
Re(nii) = ((nj|) + (n^))/2 gives a bi-fermionic distribution. The non-vanishing and 
purely imaginary quantity (n^) — (n^) describes the magnetization of the system as 
a response to the fictitious imaginary magnetic field, which is also not an observable 
quantity. 

The mean-field semi-fermionic theory of superconductivity in LNU-Hubbard 
model was developed in 14 < 37 . The basic result indicates on Bose-condensation type 
of the superconductivity with critical temperature 



2\t* 



v(2i/ - 1)(1 - vf 



1/2 

(88) 



ln(i//(l - v) 

with t* = zt. the critical temperature varies as T c ~ (i/lnT 1 i/) 1 ^ 2 for v — > and 
T c <~ \t*\/2^/2z for v — ► 1/2. We address the reader to original papers 14 ^ 37 for 
details of loop corrections and excitation spectra in the LNU-Hubbard model. 
The t-J Hamiltonian has the form 

H = -t (1 - A,-<,"i,-M,**iA'L - a],-° a i.-°) + jJ2SiSj (89) 

and can be obtained from the Hubbard model near half-filling point under condition 
t <C U. The spin exchange coupling J = At 2 jU . 

The t-J model in the form (89) does not posses the most important condition 
for the application of semi-fermionic representation. Namely, it satisfies neither the 
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particle-hole symmetry manifested for SU (2) spin systems, nor global particle-hole 
symmetry of LNU-Hubbard model. In the paper 15 Gross and Johnson proposed 
to consider generalized version of t-J model when kinetic energy is chosen to be 
particle-hole symmetric 



-* E 



[(! - a i,-a a i -vWi^jA 1 - a \,-<T a 3-o) + n i-<r a la a i,* n 3-v ( 90 ) 



where n i (y — at 0j )O .. It has been shown in 15 that the thermodynamical properties 
of original and the generalized t-J model are identical. 

The partition function of t-J model is mapped to those of generalized t-J model 
through projection operator 15 . Unlike particle-hole symmetrical spin case, where 
both empty and double occupied states have to be excluded in calculation of traces, 
only doublons (double occupied states) are unwanted in t-J model. The exclusion 
of double occupied states is not done by Schwinger-fermion procedure when the 
constraint + mi = 1 is enforced by an additional U(l) fluctuating field A 



f 

Jo 



2tt 

dA(cxp (-/3H + i\i[n^ + n a - 1])) 



but by discrete complex Lagrange multipliers (chemical potential). The "complex" 
chemical potential /i tot consists of the real part a which takes care about finite hole 
concentration and imaginary part which represents the constraint. The auxiliary 
variable a — Re/j, tot is related to the chemical potential of holes Hh through equation 

2sim>/3) =exp(/i h /3) (91) 

The half-filled case is described by the limit a — * 0. It corresponds to the Heisen- 
berg model for which semi-fermions are characterized by purely imaginary chemical 
potential /j, = inT/2. Thus, the scmi-fermionic representation is generalized for 
non-particle-hole symmetric cases as well. Another example without particle-hole 
symmetry (dynamical symmetries) will be considered in the Section devoted to 
application of semi-fermions in mesoscopic physics. 



2.4. Kondo lattices: competition between magnetic and Kondo 
correlations 

The problem of competition between Ruderman-Kittel-Kasuya-Yosida (RKKY) 
magnetic exchange and Kondo correlations is one of the most interesting problem of 
the heavy fermion physics. The recent experiments unambiguously show, that such 
a competition is responsible for many unusual properties of the integer valent heavy 
fermion compounds e.g. quantum critical behavior, unusual antifcrromagnetism and 
superconductivity (see references in 20 ). We address the reader to the review 38 for 
details of complex physics of Kondo effect in heavy fermion compounds. In this 
section we discuss the influence of Kondo effect on the competition between local 
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(magnetic, spin glass) and non-local (RVB) correlations. The Ginzburg-Landau the- 
ory for nearly antiferromagnetic Kondo lattices has been constructed in 20 using the 
semi-fermion approach. We discuss the key results of this theory. 
The Hamiltonian of the Kondo lattice (KL) model is given by 



(92) 



Here the local electron and spin density operators for conduction electrons at site 
j are defined as 

jcr a 

where f are the Pauli matrices and Cj a — ~^2 k exp(ifcj). The spin glass (SG) 
freezing is possible if an additional quenched randomness of the inter-site exchange 
Ijl between the localized spins arises. This disorder is described by 

H' ^^(SjS,). (94) 
ji 

We start with a perfect Kondo lattice. The spin correlations in KL are characterized 
by two energy scales, i.e., / <~ J 2 /sf, and ~ £f exp(— ep/ J) (the inter-site 
indirect exchange of the RKKY type and the Kondo binding energy, respectively). 
At high enough temperature, the localized spins are weakly coupled with the elec- 
tron Fermi sea having the Fermi energy £p, so that the magnetic response of a 
rare-earth sublattice of KL is of paramagnetic Curie- Weiss type. With decreasing 
temperature either a crossover to a strong-coupling Kondo singlet regime occurs at 
T ~ Ak or the phase transition to an AFM state occurs at T = Tjy ~ zl where z 
is a coordination number in KL. If TV ~ the interference between two trends 
results in the decrease of both characteristic temperatures or in suppressing one of 
them. The mechanism of suppression is based on the screening effect due to Kondo 
interaction. As we will show, the Kondo correlations screen the local order param- 
eter, but leave nonlocal correlations intact. The mechanism of Kondo screening for 
single-impurity Kondo problem is illustrated on Fig. 6 As a result, the magnetization 
of local impurity in the presence of Kondo effect is determined in terms of GF's of 
semi-fermions G(u>) by the following expression 39 : 



M(H) = S(9^ B )TJ2 (StM - SiM) = 5(. 9M B)tanh 
1 ln(ln(T/Tk)) , 



ln(T/T K ) 2\n 2 {T/T K ) 

To take into account the screening effect in the lattice model we apply the semi- 
fermionic representation of spin operators. In accordance with the general path- 
integral approach to KL's, we first integrate over fast (electron) degrees of freedom. 
The Kondo exchange interaction is decoupled by auxiliary field <j> 40 with statistics 
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Fig. 6. Kondo screening of the local moment in single-impurity Kondo problem. Dashed line 
denotes semi-fermions, solid line stands for conduction electrons. 



complementary to that of semi-fermions which prevents this field from Bose con- 
densation except at T = 0. As a result, we are left with an effective bosonic action 
describing low-energy properties of KL model at high T > Tk temperatures. 

• Kondo screening of the Neel order 

To analyze the influence of Kondo screening on formation of AFM order, we adopt 
the decoupling scheme for the Heisenberg model discussed in Section II. A. Taking 
into account the classic part of Neel field, we calculate the Kondo-contribution to 
the effective action which depends on magnetic order parameter A/": 



l-n(A0 



Mq)| 2 - 



(96) 



where a polarization operator II (AO casts the form 

1 



II(A0=p(0)ln( 



3h(/W) 



1 



<>(£) 



(97) 



where p(0) is the density of states of conduction electrons at the Fermi level and 
the Kondo temperature Tk = £f cxp (— l/(p(0) J)). Minimizing the effective action 
A(4>, AO with respect to classic field AO the mean field equation for Neel transition 
is obtained (c.f. with (61)) 



JV = tanh -§ 



2T 



1 



cosh 2 (/?/ Q AA/2) 



In (T/T K ) cosh 2 (/3/ Q A0 



(98) 



As a result, Kondo corrections to the molecular field equation reduce the Neel 
temperature 



• Kondo enhancement of RVB correlations 
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Applying the similar procedure to nonlocal RVB correlations, we take into account 
the influence of Kondo effect on RVB correlations 



E 



cosh/?(7 k A) 



n(/ q A)=p(0)ln(^) 
- 1 + 7 k A tanh(/3/ k A) 



+ 



^ +q + (V2/?) 2 ' 



(99) 



Here = e(k) — ep. Minimizing the effective action with respect to A we obtain 
new self-consistent equation to determine the non-local semi-fcrmion correlator. 



tanh 



la A 



+ a s 



J q A 



T\n{T/T K ) 



(100) 



It is seen that unlike the case of local magnetic order, the Kondo scattering favors 
transition into the spin-liquid state, because the scattering means the involvement 
of the itinerant electron degrees of freedom into the spinon dynamics. 

2.5. Semi-fermionic representation for spin glass models 

In this Section we sketch the technicalities related to the application of Popov- 
Fcdotov method to disordered spin systems. We address the reader to a review 41 
on theoretical concepts and experimental facts on spin glasses. 
We consider Ising/Hciscnberg spin glass model 



H 



<ij> 



(101) 



with gaussian distributed random interaction and A = corresponds to 
Sherrington-Kirpatrick (SK) model 42 , while A = 1 corresponds to isotropic Heisen- 
berg model. 

Following the procedure described in 43,44 we use fermionic representation of 
spin operators employing the replica trick 43 : 

Oi(T) -> #(t), o=l..n. (102) 

As a result we obtain the expression for the average of n-th power of the partition 
function for isotropic model: 

(Z n ) av = J] / dlijPiln) J! D &U ex P {[ dr 



(103) 



<i,j> a=l 

Assuming P{hj) to be a gaussian distribution P(Iij) ~ exp(— If N/{2I 2 )) we in- 
tegrate over Iij. As a result, four-spin, or, eight-fermion term appears in effective 
action. To decouple the quantity 

*;W) = E E ^°o(^<^w^( T >;i^') ( 104 ) 

i aa'77' 
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the Q-matrices are introduced: 



(Z n )av = J D[Q]cxp ^((3I) 2 NTr[Q 2 } + 



(105) 



and (Gq )aa' = 5aa' {iu n — oigh/2). For Ising model only X^ b should be considered. 

The next step is to consider the second decoupling procedure for lr[QX). For 
this sake, following 43 ' 44 we consider Q as a constant saddle-point matrix. 

QTp = Q, Qsf = Q (106) 

The replica-global term S 1 ") 2 is decoupled by a replica independent Hubbard- 
Stratonovich field, while the replica local product (S£) 2 requires a replica-dependent 
decoupling field . Although it is well known that both SK and isotropic vector 
model of infinite range spin glass are unstable with respect to Parisi's broken sym- 
metry permutation symmetry 45 , we shall consider first the replica symmetrical 
solutions existing only near the phase transition point. 

The last step is to integrate over all Popov-Fedotov fermions taking the limit 
n — > and calculate free energy per site 

f = T ^( 1 - <Zn> ^ < 107 > 

and find the extremum of / with respect to q and q. 

We emphasize, that in the framework of quantum theory treatment we do not 
have a freedom to choose q freely. In general, the equations for q and q are the set 
of coupled equations. 

a) Sherrington-Kirkpatrick model. 

After applying the scheme described above to Ising spin-glass model (A = 0) 
and some straightforward calculations one gets the following expressions: 

exp Q(/3/) 2 Tr[QA]) = £ ]\ £ exp f £ (*, !/>£»] (108) 



a, a 

G p+oo 

2 



(-x 2 /2)... 



/>u r-\-oc 

where H(z 1 y a ) = Iy/qz + lyjq — qy a and / ... = 1/V2k / (ixexp 

J x J — oo 

Then, the integration over ip a becomes gaussian and can be performed explicitly. 

Performing the integration over y a the following expression for the free energy 
is obtained: 



\2 r G 



Pf = ((! - «) 2 - (! " if) - I M 2 cosh(/3/^)] (109) 
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resulting in the equations: 



q= / tanh 2 (/?/V^), 3 = 1- 

J z 

Thus, the SK results are exactly reproduced by Popov-Fedotov method. 



(110) 



b) Isotropic vector spin-glasses. 

Applying the same procedure to isotropic vector spin-glass model, we obtain 
the set of equations for diagonal q and non-diagonal q saddle point elements of Q 
matrix 43 < 44 : 



Q = 



2 



f 

Jo 



drexp — 



2q(piy 



b 2 + tanh(r) [r-b 2 /r] 



q = 



3 g 3 / 2 (/?/) 3 v / 2Wo 



dr exp 



2q(/3I)< 



l + b 2 tanh(r)/r 
r 2 + rtanh(r)[& + 2] 



l + b 2 tanh(r)/r 



(111) 



(112) 



with b = \J q — q(3I. These equations are the generalization of equations Eq.110 for 
replica-symmetric solution of isotropic Heisenberg spin glass model. We address the 
readers to original papers 43 ' 44 for discussion of stability (de Almeida-Thouless) line 
of the vector model and broken replica solutions in the framework of semi-fermionic 
approach. 



2.6. Kondo lattices with quenched disorder 

In this Section we consider the Kondo Lattice model (92) with random RKKY 
interaction. The randomness is associated with e.g. the presence of non-magnetic 
impurities in magnetic compound resulting in appearance of random phase in the 
RKKY indirect exchange. In this case the spin glass phase should be considered. As 
it has been shown in 16 and 20 , the influence of static disorder on Kondo effect in 
models with Ising exchange on fully connected lattices (SK model) can be taken into 
account by the mapping KL model with quenched disorder onto the single impurity 
Kondo model in random (depending on replicas) magnetic field. It allows for the 
self-consistent determination of the Edwards- Anderson qsA order parameter given 
by the following set of self-consistent equations 



Ih(T/Tk) W{T/T K )) ' 

q = L tanh \YTwrw^wm) + wmr) ) ■ (113) 

Here q = qEA and q are nondiagonal and diagonal elements of Parisi matrix re- 
spectively. Therefore, the Kondo-scattering results in the depression of the freezing 
temperature due to the screening effects in the same way as the magnetic mo- 
ments and the one-site susceptibility arc screened in the single-impurity Kondo 
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problem when Ising and Kondo interactions are of the same order of magnitude. 
Let's now briefly discuss the fluctuation effects in Kondo lattices. The natural way 




Fig. 7. Fcynman diagrams for nonlocal excitations associated with the overlap of Kondo clouds. 
Wavy line denotes Kondo "scmi-Bosonic" field, solid line stands for fermions, dashed line denotes 
semi-fermions. 

to construct the fluctuation theory is to consider the non-local dynamical Kondo 
correlations described by the field 0(q, u>) (see Fig. 7). In fact, the non-locality of the 
"semi-Bosonic" field is associated with an overlap of Kondo clouds 20 and responsi- 
ble for a crossover from the localized magnetism to the itinerant-like fluctuational 
spin-liquid magnetism. The temperature dependence of static magnetic susceptibil- 
ity becomes nonuniversal in spite of the fact that we are in a region of critical AFM 
fluctuations which is consistent with recent experimental observations (see 20 for 
details). 

3. Application of semi-fermions in mesoscopic and nano-physics 
3.1. Kondo effect in quantum dots 

The single electron tunneling through the quantum dot 46 has been studied in great 
details during the recent decade. Among many interesting phenomena behind the 
unusual transport properties of mesoscopic systems, the Kondo effect in quantum 
dots, recently observed experimentally, continues to attract an attention both of 
experimental and theoretical communities. The modern nanoscience technologies 
allow one to produce the highly controllable systems based on quantum dot devices 
and possessing many of properties of strongly correlated electron systems. The quan- 
tum dot in a semiconductor planar heterostructure is a confined few-electron system 
(see Fig. 8) contacted by sheets of two-dimensional gas (leads). Junctions between 
dot and leads produce the exchange interaction between the spins of the dot and 
spins of itinerant 2D electron gas. Measuring the dc I — V characteristics, one can 
investigate the Kondo effect in quantum dots under various conditions. 

Various realizations of Kondo effect in quantum dots were proposed both theo- 
retically and experimentally in recent publications (see e.g. 47,48 for review). In order 
to illustrate the application of scmi-fermionic approach we discuss briefly electric 
field induced Kondo tunneling in double quantum dot (DQD). As was noticed in 
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49 , quantum dots with even AT possess the dynamical symmetry SO (4) of spin ro- 
tator in the Kondo tunneling regime, provided the low-energy part of its spectrum 
is formed by a singlet-triplet (ST) pair, and all other excitations are separated from 
the ST manifold by a gap noticeably exceeding the tunneling rate T. A DQD with 
even J\f in a side-bound configuration where two wells are coupled by the tunneling 
v and only one of them (say, I) is coupled to metallic leads (L,R) is a simplest 
system satisfying this condition 49 > 50 . Such system was realized experimentally in 
Ref. 51 . As it was shown in 23 the Shrieffer- Wolff (SW) transformation, when applied 



Left lead (L) 



W 




w 



Right lead (R) 

a) 



K 



b) 



Fig. 8. (a) Double quantum dot in a side-bound configuration (b) co-tunneling processes in biased 
DQD responsible for the resonance Kondo tunneling. 



to a spin rotator results in the following effective spin Hamiltonian 

Hint = J a oL'.ft.fs c \ aa c k'a'o 
kk' ,aa' —L^R 

+ E (jL'Si A '+J^RiA')r^4 ac7 c kWal f{f AI (114) 

kk' ,aa' AA' 

where the c-operators describe the electrons in the leads and /-operators stand for 
the electrons in the dot. The matrices S and R (d=x,y,z) are 4x4 matrices 
defined by relations (A. 13) (see Appendix A) and J s — J ss , J T = J TT and J ST 
are singlet, triplet and singlet-triplet coupling SW constants, respectively. 

Applying the semi-fermionic representation of 50(4) group introduced in Section 
I, we started with perturbation theory results analyzing the most divergent Feynman 
diagrams (Fig. 9) for spin-rotator model 23 . 

Following the "poor man's scaling" approach we derive the system of coupled 
renormalization group equations (see Fig. 9) for effective couplings responsible for 
the finite bias eV transport through DQD. 

As a result, the differential conductance G(eV, T)/Gq ~ \ Jlr\ 2 IS shown to be 
the universal function of two parameters T/Tk and eV/Tx, Go = e 2 /nh: 

G/G - In" 2 (max[(eV r - 6), T]/T K ) (115) 

where Tk is a non-equilibrium Kondo temperature 23 . Thus, the tunneling through 
singlet DQDs with the singlet /triplet splitting 5 = Et — Es 3> Tk exhibits a peak 
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Fig. 9. Right panel: leading (b,d) and next to leading (c,e) parquet diagrams determining renor- 
malization of J s (a). Solid lines denote electrons in the leads. Dashed lines stand for semi-fermions 
representing double electron states in the dot. Left panel: irreducible diagrams contributing to RG 
equations. Hatched boxes and circles stand for triplet-triplet and singlet-triplet vertices respec- 
tively. 



in differential conductance at eV w 6 instead of the usual zero bias Kondo anomaly 
which arises in the opposite limit, 8 < T K . Therefore, in this case the Kondo effect 
in DQD is induced by a strong external bias. 

The scaling equations can also be derived in Schwinger-Keldysh formalism (see 
12 and also 19 ) by applying the "poor man's scaling" approach directly to the dot 
conductance. Thus, applying the semi-fermionic representation to the local (zero- 
dimensional) problem with dynamical symmetry 48 we developed regular perturba- 
tion theory approach and obtained the scaling equations. 

3.2. From double quantum dot to spin chains and ladders 

In this Section we show the application of semi-fermionic representation to ID sys- 
tems (spin chains and ladders). The semi-fermions give a powerful tool for fermion- 
ization of the spin system which can be bosonized as a next step in ID. The bosoniza- 
tion procedure performed by means of Hubbard-Stratonovich transformation 52 ' 53 ' 54 
is known in the literature on Tomonaga-Luttingcr model as "functional bosoniza- 
tion" (see also 55 > 56 ). It is known to be asymptotically exact in the long-wave limit. 
The semi-fermionic representation allows us to apply the functional bosonization 
methods to ID spin systems 57 . 

Let us consider, as an example Spin Ladder (SL) and Spin Rotator Chain (SRC) 
models, Fig. 10 58 . Generic Hamiltonian for spin systems under consideration is the 
Heisenberg-type spin 1/2 ladder Hamiltonian 

H^ = J t Sil ■ So + Jj ^ S -' S J« ( 116 ) 

(il,»2) a (ia.ja) 

Here index a = 1, 2 enumerates the legs of the ladder, and the sites (zl, i2) belong 
to the same rung (Fig. 10a). 

A chain of dimers of localized spins illustrated by Fig. 10b is described by the 
simplified version of this Hamiltonian 

H SRC = J t J2 S<i • S j2 + Ji Y Bil ' Sj! (117) 

(n,i2) (ij) 
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The geometry of alternate rungs is chosen in a system (Fig. 10b) to avoid exchange 
interaction between spins Si2 and Sj%. 

The transverse coupling may emerge either from direct exchange (in case of local- 
ized spins) or from indirect Anderson-type exchange induced by tunneling (similarly 
to the case of quantum dots). In the latter case the sign of Jt is antiferromagnetic, 
in the former case it may be ferromagnetic as well. The same is valid for <//. 



(a) (b) 
Fig. 10. Spin Ladder (a) and Spin Rotator Chain (b) 



As it is shown in 2 <> 58 >° s > SO (4) is an intrinsic property of S=l/2 spin ladders 
and decorated spin chains. 

To describe the elementary excitations in SRC we use the semi-fermionic repre- 
sentation for SO (4) group (see Section I). Then the single-site Hamiltonian may be 
represented in a form 

Hi = -EsTflhs + Kflfa - fl-Ja) (H8) 

where h is uniform magnetic field applied in z direction and Est is a singlet / triplet 
splitting. 

Let us consider the XXZ-SRC model 27 . The Hamiltonian of this model is written 
in terms of the generators of SO (4) group (see Appendix B) as follows 

H = i Ji E ( s t s i + s t R j + s ^ R i + A i s * s 3 + 2S ! RZ j}) ■ ( n9 ) 

(ij) 

The S-S part of this Hamiltonian describe the 5* = 1 chain, with the Haldane 
gap in the excitation spectrum (see, e.g., 60 ' 59 ). The question is, how do the S-R 
interaction modifies the gap. We consider the case of FM dimers, when the triplet 
is the ground state. In this case one has one more gap mode, where the gap equals 
Jt- This mode is coupled to Haldane branch only via S-R exchange terms in (119). 

The spin liquid fcrmionization approach adopted here is a convenient tool for 
description of Haldane spectrum. Unlike the 5=1/2 model, where the spin-liquid 
state is easily described by globally U(l) invariant modes TijTji = J2 a \fiafj<?\ 2 ' m 
case of S — 1, one deals with variables which effectively break this symmetry. One 
can rewrite the effective Hamiltonian of SRC model with A = in a form 

H =\jiJ2 (Afn + fllffl) BfBf + fkflcfBf + Bfcff^fa, (120) 

where B° s = f 0j + fsj , C° 5 = f j - fsj- The terms in the first line of Eq. 
(120) describe coherent propagation of spin fermions accompanied by a back-flow 
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on neutral fcrmions, whereas the terms in the second line are "anomalous" (they 
do not conserve spin fermion number). For example the propagator (S^SJ) con- 
tains anomalous components //i/jj/jo/io - » F*. 1I F jii0 o along with normal ones 
/ii/ji/jo/io- Here F^.aa' = /ja/»A'- The first term in (120) describes the kinetic 
energy spinon excitations, and two last anomalous term breaking U(l) symmetry 
are responsible for the Haldanc gap. 




Fig. 11. Lowest order contributions to anomalous propagator. 



To reveal the contribution of dynamical symmetry on the Haldane gap, one have 
to note that these terms enter both in countcrflow in the first term and in gauge 
symmetry breaking terms in the second line. In spin 1 ladder the countcrflow term 
~ fiofjo predetermines the width of spinon band described by the first line of Eq. 
(120). Apparently, the one more channel (tripet/singlet transitions in B os ) enhances 
this effect, because in this case the local constraint imposes more restrictions of 
phase fluctuations. 

The gap itself is due to anomalous correlations described by the second line of 
Eq. (120). Here the appearance of second channel of spinless excitations results in 
formation of even and odd operators B® and Cj . The Haldane gap closes when 
the |0) and \S) states are degenerate (the odd operator C® s nullifies the anomalous 
terms responsible for its formation). This means that appearance of OS channel 
favors closing of the Haldane gap. 

In a strong coupling case of J t » Ji both above trends may be considered at 
least in the lowest order of a perturbation theory. In case of spin ladders 61 ' 62 the 
1-st and 2-nd-order in g = Ji/J t anomalous diagrams are represented in Fig. 11. 

4. Epilogue and perspectives 

In this paper, we demonstrated several examples of the applications of semi- 
fermionic representation to various problems of condensed matter physics. The list 
of these applications is not exhaustive. We did not discuss, e.g., the interesting de- 
velopment of SF approach for 2D Ising model in transverse magnetic field, functional 
bosonization approach based on semi-fcrmionic representation applied to S = 1/2 
chains, application of SF formalism to mesoscopic physics 24 etc. Nevertheless, we 
would like to point out some problems of strongly correlated physics where the 
application of SF representation might be a promising alternative to existing field- 
theoretical methods. 



3, 2008 10:51 WSPC/INSTRUCTION FILE ijmpb'kiselev 



Semi-fermionic representation for spin systems ... 33 

Heavy Fermions 

• Crossover from localized to itinerant magnetism in Kondo lattices 

• Quantum critical phenomena associated with competition between local 
and nonlocal correlations 

• Noncquilibrium spin liquids 

• Effects of spin impurities and defects in spin liquids 

• Crystalline Electric Field excitations in spin liquids 

• Dynamic theory of screening effects in Kondo spin glasses. 

Mesoscopic systems 

• Noncquilibrium Kondo effect in Quantum Dots 

• Two-channel Kondo in complex multiple dots 

• Spin chains, rings and ladders 

• Nonequilibrium spin transport in wires 

Summarizing, we constructed a general concept of semi-fermionic representation for 
SU (N) and SO(N) groups. The main advantage of this representation in application 
to the strongly correlated systems in comparison with another methods is that the 
local constraint is taken into account exactly and the usual Feynman diagrammatic 
codex is applicable. The method proposed allows us to treat spins on the same 
footing as Fermi and Bose systems. The semi-fermionic approach can be helpful for 
the description of the quantum systems in the vicinity of a quantum phase transition 
point and for the nonequilibrium spin systems. 
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Appendix A. SO(N) group 

Lie algebras o(n) are defined on a basis of |n(n — 1) operators of infinitesimal 
rotation 

d d 

D a p = -Dp a = xp — Xa Q^> {a < (3 = l,2,...n) (A.l) 

which possess the following commutation relations 

[DapD^v] = {5 a ^Dp v - 5 av Dp^ - Sp^D^ + Sp^Da^). (A. 2) 

(see, e.g., 63 < 64 ). 
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We consider o(4) algebra as a first example of dynamical symmetries 63 . The 
antisymmetric tensor D a p for o(4) algebra acts in 4-dimcnsional space. It can be 
parameterized in terms of two vectors L and M as follows 



/0L 3 -L 2 MA 
Li M 2 
M 3 

V o / 



(A.3) 



where the infinitesimal operators of 5*0(4) group 64 in 4-dimensional space (x, y, z, t) 
are given by 

( d d\ T ( d d\ T ( d d 
L\=i \ y-z z— , L 2 = i [ z- x— , L 3 = i[ x- y— 

\ Oz ay ) \ ox oz J \ oy ox 

^^m-**)* M2 = l {%- y di)> M3= T^-^> (A - 4) 

From here we deduce the commutation relations: 

[Lj,L k ] = iejkiLi, [Mj,M k ] = ie jk iL u [Mj,L k ] = ie jki Mi. (A.5) 
The two Casimir operators consist of nontrivial 

L 2 + M 2 = 3 (A.6) 
and trivial orthogonality condition 

L • M = M • L = 0. (A.7) 

It is known, that three generators of SU(2) group together with the Casimir 
operator L 2 define a sphere 5*2 where each state is parameterized by two angles. 
The coherent states of SU (2) group may be constructed 65 by making a standard 
stereographical projection of the sphere from its south pole to the complex plane 
z. The space of generators of SO(4) group is 6- dimensional, while 2 constraints 
determine 4-dimensional surface, where each state is characterized by four angles. 
The stereographical projection of this surface on a 4-D complex hypcrplanc allows 
to construct coherent states for SO (4) group. 

The commutation relations (A.5) can be transformed into another form by mak- 
ing the linear transformation to the basis 

J, = i!±Mi, K, = i^ (A.8, 

giving more simple commutation relations 

[Jj, Jfc] = iejkiJi, [Kj,K k ]=ie jk iKi, [Kj,J k ]=0. (A. 9) 

The operators Li, Mi as well as Ji, Ki form the elements of the Lie algebra o(4). The 
operators (Ji, J2, J3) and {K\, K 2 , if 3) are separately closed under commutations, 
each describing a subalgebra of o(4), namely o(3) = u{2). The Lie algebra o(4) 
is the direct sum of two o(3) algebras. This splitting of the o(4) algebra into two 
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o(3) subalgebras is directly associated with the local isomorphism between the Lie 
group SO(4) with the direct product group SU(2) x SU(2). The triads (Ji, J 2 , J 3 ) 
and (Ki, K 2 , K 3 ) each form proper ideals 64 in o(4), and the Lie algebra o(4) is 
semi-simple. 

The symmetry group of spin rotator is defined in close analogy with the above 
construction, but all rotations are performed in a spin space. The triplet/singlet pair 
is formed in a simplest case by two electrons represented by their spins s — 1/2. Let 
us denote them as si and S2 . The components of these vectors obey the commutation 
relations 



[sij, sifej = lejkisu, 



[S2j, S 2 k] = iejklS2l, 



[Slj, S 2 k\ 







(A.10) 



In similarity with (A. 9) these vectors may be qualified as generators of o(4) algebra, 
which represents a spin rotator. Then, the linear combinations 



Sli 



S2i 



Ri 



Sli 



S2i 



(A.11) 



arc introduced in analogy with (A. 8), which define 6 generators of SO(4) group 
possessing the commutation relations (A. 5). These generators are represented in 
terms of the Pauli-like matrices as follows 



S+ = V2 



/0 1 0\ 
10 


\o o o o/ 



V2 



fO 0\ 
10 
10 

Vo o o 0/ 



s z = 



R+ = V2 



where the ladder operators S ± 
(A. 7) now acquire the form 



no 





o\ 














-1 





\0 





0/ 


/0 





!\ 




















\0 


-1 


0/ 



R z = - 



/0 0\ 
' 00 1 


Vo i o oj 



R~ = V2 



/0 \ 

-1 

\1 / 



(A.12) 



(A.13) 



S x ±iS y , R ± = R x ±iRy. The constraints (A.6), 



S • R = 0, S 2 + R 2 



3. 



By construction S is the operator of the total spin of pair, which can take values 
S = for singlet and S — 1 for triplet states. The second operator R is responsible 
for transition between singlet and triplet states. Thus we come to the dynamical 
group SO (4) for spin rotator. 
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Similar procedure is used for the 50(5) group. The corresponding o(5) algebra 
has 10 generators D a p = —Dp a (A.l) satisfying commutation relations (A. 2). These 
10 generators may be identified as 3 vectors and a scalar in a following fashion 

-5" R x P x \ 

gx R y py 

R z P z 
A 
J 



/o s z 





—I 



\ 



(A.14) 



where the operators S, R, P and the scalar operator A obey the following commu- 
tation relations 



[Sj,Sk] — iejkiSi, [Rj,Rk] = ie jk iSi, 
[R 3 ,Sk] = iejkiRi, 
[Pj,A]=iRi, 



[Pj,S k ] = iejkiPi, 



[A,R 3 ]^tP„ 



[Pj,Pk] = iejkiSi, 
[Rj,Pk] = i5 jk T, 
[A,Sj} = 0. (A.15) 

The operators R and P are orthogonal to S, while the Casimir operator is 
K, = S 2 + R 2 + P 2 + A 2 =4. These operators act in 10-D spin space, and the 
kinematical restrictions reduce this dimension to 7. Similarly to 50(4) group, the 
vector operators describe spin S=l and transitions between spin triplet and two 
singlet components of the multiplet, whereas the scalar A stands for transitions 
between two singlet states. 

The group 50(5) is non-compact, and the parametrization (A.14) is not unique. 
As an alternative, one may refer to another representation of D a p used, the theory 
of high-T c superconductivity 66 : 

/ \ 


-S z (A.16) 

S v -S x 

- 71 x ) -i(7T+ - TT y ) — i(7T+ - IT z ) 0/ 

with 10 generators identified as a scalar Q and three vectors 5, ir and 7?t standing 
for the total charge, spin and it triplet 5 = 1 superconductor order parameter. Both 
representations (A.14) and (A.16) are connected by the unitary transformation. 



TTy + 7T, 
TTz+TT 

Q 



Appendix B. Dynamical Symmetries in Spin Rotator Chain Model 

The Spin Ladder and Spin Rotator Chain Hamiltonians are given by 

(tl,t2) a (ia,ja) 

H SRC = J t sn-sa + Ji^sa-Sj! (B.2) 

(n,i2) (ij) 

We start with diagonalization of the Hamiltonian of perpendicularly aligned 
dimer. The 50(4) symmetry stems from the obvious fact that the spin spectrum 
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of a dimer {il, 12} is formed by the same singlet-triplet (ST) pair as the spin spec- 
trum of DQD studied in the previous section. This analogy prompts us a canonical 
transformation connecting two pairs of spin vectors, {sji,s,2} and {Sj,Ri}: Two 
sets of spin operators are connected by a canonical transformation 

Si + R,i Si — R,i . . 

s*i = ^ ' Si2 = 2 ' ^ ' 

Then the Hamiltonian Hi of a single dimer i is the same as the Hamiltonian (114) 
of DQD. The total spin of a dimer is not conserved in a spin chain, so the dynamical 
symmetry of individual rung is revealed by the modes propagating along the chain 
62 . Applying the rotation operation to the Hamiltonians (116), (B.2), we transform 
them to a form 



Here 



H = Ho+ H int (B.4) 



Ho — y^^Hj, 



is common for both models. It is useful to include the Zeeman term in Hi, 

Hi = \ (E S B% + E T Sj) + hS lz . (B.5) 

We confine ourselves by a charge sector iVj = 2 and omit the Coulomb blockade 
term for the sake of brevity. The interaction part of SL Hamiltonian transforms 
under rotation (B.3) to the following expression 

The interaction part of the SRC Hamiltonian is 

KtnT = \ji E( S *^ + + R ^) (B-7) 

Now we see that both effective Hamiltonians belong to the same family. In all 
three case one may transform initial ladder or "semi-ladder" Hamiltonians into 
really one-dimensional spin-chain Hamiltonians, which, however, takes into account 
the hidden symmetry of a dimer. The effective Hamiltonians (B.6), (B.7)) contain 
operators R describing the dynamical symmetry of dimers. It is clear that this 
dynamical symmetry make the spectrum of this Hamiltonians more rich than that 
of a standard Hciscnbcrg chain. Like in many other cases, rotation transformation 
eliminates the antisymmetric combination of two generators. 

The transformation (B.3) reveals the hidden symmetry of well known spin 1/2 
ladder (B.6). It maps the ladder hamiltonian onto a pair of coupled chain Hamiltoni- 
ans: one is conventional spin 1 Haldanc chain , another is unconventional pseudospin 
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1 Haldane chain. Spin and pseudospin are coupled dynamically by the commutation 
relations and kinematically by the local Casimir constraint 

S 4 2 + R2 = 3. (B.8) 

One may also compare the Hamiltonian (B.6), (B.7) with effective Hamiltonian of 
spin-1 chain, which arises after decomposition of spin-one operators into pair of 
spin 1/2 operators, Sj = Sj + r^. This decomposition operation transforms initial 
Hamiltonian into a form similar to H c but for spin-one- half operators Sj, r^. The 
difference between two cases is that these effective spins commute unlike operators 

Hi. In other words, the difference is that the local symmetry of spin-one chain 
is 50(3) whereas the local symmetry of SRC is 50(4). The spin rotator chain (B.7) 
is in some sense intermediate between spin chains and spin ladders. In this cases 
the spin-pseudospin symmetry is obviously broken by the cross terms 2SiR.j, and 
our aim is to find the specific features of this novel model systems. 

In all cases the simplified versions of Heisenbcrg Hamiltonians may be consid- 
ered. The simplified SL models are well known 61 . The simplified anisotropic versions 
of the Hamiltonian (B.7) have the following forms: 

Ising-like SRC model: 

H = \ Jl E W + 2SZ i RZ j + R t R j) (B-9) 

to) 

Anisotropic SRC model: 
H = \ J iY, [^ S 7 + S t R J + S i R t + R t R j) + + 2S?R* + RfRj)] 

m 

SRC in strong magnetic field: SO (4) group is reduced to SU(2) group in strong 
magnetic field, when the Zeeman splitting exactly compensates the exchange gap 
in a single dimer, ho = \Et — E$\. Then at low T, the states |i0) and \i — 1) are 
quenched, and only two components, R^ survive in the whole manifold . As a result, 
the Hamiltonian (B.7) is mapped onto a XF-modcl for spin one half: 

H xy = \ji E(^ +i? 7 + H - c -)- (B.10) 

This means that starting from a singlet ground state for J t = E T — E s > 0, one 
may induce development of spin liquid-like excitations by applying strong magnetic 
field. In a near vicinity of this point of degeneracy, H int acquire the feature of XY 
model in transverse magnetic field. 
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